Generation of NPT Entanglement from Nonclassical Photon Statistics 
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With a product state of the form pin = pa® |0)66{0| as input, the output two-mode state pout of 
the beam splitter is shown to be NPT whenever the photon number distribution (PND) statistics 
{p(na)} associated with the possibly mixed state pa of the a-mode is antibunched or otherwise 
nonclassical, i.e., if {pijia) } fails to respect any one of an infinite sequence of classicality conditions. 
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Considerable progress has been achieved in recent 
years towards an understanding of entanglement in the 
context of continuous variable bipartite systemsQ, 0, 0, 
0]. Phase space descriptions and Gaussian states have 
played a significant role in this progress. 

While the quantum uncertainty principle places re- 
strictions on the moments of the phase space variables, 
separability of a state places additional demands on these 
moments. Gaussian states are fully determined by their 
first and second moments, or variances, and so also is 
the issue of their separability. But non-Gaussian states 
which satisfy the separability demand on their variances 
could still be entangled by virtue of violation of one of 
the separability requirements on their higher order mo- 
ments. It is only relatively recently that interest in this 
direction has started emergingj^. 

There is a useful connection between nonclassicality 
and inseparability, and the beam splitter plays an im- 
portant role as a bridge between these two attributes. 
Asboth et alJ^ have shown that the output of a beam 
splitter whose input is a product state pin = Pa® |0)bb(0| 
is entangled if and only if the single-mode state pa at the 
input is nonclassical. [Thus any measure of entanglement 
of the output state in this configuration is a computable 
measure of nonclassicality of pa, the entanglement po- 
tential {EP)0]. For use in quantum information tasks, 
however, it is important to know if the entanglement gen- 
erated in this manner is NPT or PPT. 

Two canonical manifestations of nonclassicality have 
been extensively studied in the quantum optics context: 
(1) squeezing^, and (2) antibunching or sub-Poissonian 
fluctuation's', which is a particular manifestation of non- 
classical photon statistics; these are respectively non- 
classicalities of the phase-sensitive and phase-insensitive 
types. Whereas the former has been well explored as a 
source of entanglement in the context of Gaussian states, 
the same cannot be said in respect of the latter. 

The principal aim of this Letter is to show that with 
Pin = Pa ® |0)(,6(0| as input, the two-mode state after 
the beam splitter is definitely NPT if pa is nonclassi- 
cal of the phase-insensitive type, that is if the photon 
number distribution (PND) statistics {p{na) } associated 
with pa of the a-mode is antibunched, or possesses any 
other higher order nonclassicality. We prove this result in 



two stages: first for the restricted case of antibunched in- 
put Pa wherein we exhibit a surprisingly simple witness; 
then in the case of arbitrary PND we take advantage of 
a complete characterization of nonclassicality based on 
the theory of moments^, result and mode of 

proof extends the domain of effectiveness of the partial 
transpose criterion for separability considerably beyond 
the traditional Gaussian regime. 

Consider therefore a bipartite system, each part being 
a single mode radiation field, with respective creation 
and annihilation operator pairs a\ a and 6^, b acting on 
Hilbert spaces 7i", Ti,^ . Their only nonvanishing commu- 
tators are [a, a^] = [b,b''] = 1. The Fock or photon num- 
ber states for the two modes provide convenient ONB's 
for W, ?{'' respectively: 

K) = K!)"'/'(a^)"lO)<,, 

\nb) = (nb!)-'/'(6t)">),, na,n, = 0,1,2,... (1) 

The products \na,nii) = \na) (E" form an ONE for 
(g) H^. The partial transpose of a bipartite state p is 
the operator p on Ti,^ defined in this ONB thus: 



{na,nb\p\na ,nb) ^ {na,nf,\p\n^,nb) . 



(2) 



We have at our disposal this key fact 11]: if the partial 
transpose p is not a valid density operator (NPT), i.e., if 
p ^ 0, then p is definitely an entangled state. 

Clearly, a way to show that p is NPT is to exhibit an 
operator A of the bipartite system such that the 'expec- 
tation value' of the positive operator A in p is negative. 
We can then conclude: 

Tr (pabA'iA) < 0, for some i ^ p is NPT. (3) 

Alternatively, and in a sense more directly, we may be 
able to find some principal submatrix of the matrix 
{na,nb\p\na ,nb ) representing p, i.e., a submatrix formed 
by intersections of any subset of rows of this matrix and 
the corresponding columns^ such that this submatrix is 
not positive definite. Then again we can conclude: 

A principal submatrix of p ^ p is NPT. (4) 

We use both strategies (3) and (4) in what follows. 
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Nonclassicality and Photon Statistics: For a moment let 
us restrict attention to states pa of the a-mode alone. The 
separation of its (quantum) states into so called 'classical' 
and 'nonclassical' types is based on the diagonal coherent 
state representation|l2j of pa'- 



Pa 



n d Za (t>{Za)\Za){Za\ 



E 



The weight function (t>{za) is in general a singular distri- 
bution of a certain class. While it is real on account of 
hermiticity of pa and normalized to unity, it is in gen- 
eral not pointwisc nonnegative over the complex plane C. 
Classical states are identified as follows. 



Pa 'classical' <^ 0(za) > for all e C. 



(5) 



All other states are declared 'nonclassical'. 

All information about any state pa is captured by 
the expectation values of all possible hermitian observ- 
ables A of the system, namely Ti{paA). While any 
operator A can be written in the normal ordered form 

^ = Sjfe=o '^ife^^"'^'^' '^^ limit ourselves to observablcs 
that are 'phase invariant' i.e., remain unchanged under 



for all a. These are of the 



the map d e* 

simple form A = '^J^QCjCi^^d^ , involving a single sum, 
and are just functions of the number operator Na = a^'a. 

We can now ask for the way in which nonclassicality 
can manifest itself if only measurements of such operators 
are carried out. Clearly these expectation values involve 
only the phase-averaged information contained in (t>{za)'. 

/■oo / oo \ 

Tr(p,i) = dIaP{Ia) \Y.^CjIJ j , 

P{Ia) = H^"'^'^'') ' l^dIaP{Ia) = l. 

While P{Ia) is of course real, it may not be pointwise 
nonnegative. For a given P(Ia), the probabilities consti- 
tuting the photon number distribution (PND) are 



p{na) = {na\pa\na) 



dIaP{Ia)- (6) 



Whatever the nature of P(/a) may be, these probabilities 
are real nonnegative and add up to unity. 

We now present the following definition which will suf- 
fice for our purposes 9]: 

Pa is classical in phase invariant sense <^ P{Ia) > 0. 

(7) 

Otherwise pa is 'nonclassical' in the phase invariant sense. 
This definition is a coarse-grained version of the earlier 



definition (5) given in terms of 4>{za). Namely for the 
three mutually exclusive situations 



(i) 0(^a)>O (hence P(/a) > 0) , 
{ii) 0(^a)^O but P(/a)>0, 
(zm) P{Ia) t (hence 4){Za) t 0) : 



(8) 



while definition (5) would describe {i) alone as 'classi- 
cal' and (m), [Hi) as two levels of 'nonclassicality', the 
phase-invariant definition (7) clubs {i) and {ii) together 
as 'classical' and only {Hi) as 'nonclassical'. 

One very familiar and well known signature of phase 
invariant nonclassicality is the antibunching condition 



dIaP{Ia){Ia - { Ia)f = (A/a)' < 0, (9) 



which definitely implies P{Ia) ^ 0. But this is by no 
means the only such signature. A complete set of nec- 
essary and sufficient conditions expressing the content of 
classicality, i.e., P{Ia) > 0, in terms of PND {p{na)} will 
be described later in this Letter.. 

Transpose Operation and Expectation Values for Oper- 
ators: While transpose operation on the density opera- 
tor transcribes into 'momentum reversal' in the Wigner 
phase space description^, it is more convenient for our 
present purposes to transfer this operation directly onto 
operators and their expectation values. The transpose 
Pa of the density operator pa (of the a-mode) is defined 
in the Fock basis through (n|p^|n') = {n'\pa\n). Now the 
key point here is that in the Fock basis the operators d, d^ 
are represented by real matrices^ so the transposition op- 
eration for the matrix of d^'d™ coincides with hermitian 
conjugation, and we have Tr (p^d't-'d'^) — Tr (pad^'^d-'). 
A linear combination j, Cjjjdt^ d'^, with complex coef- 
ficients { Cj.fe }, would thus be taken under the transpose 
map to Y^ - f^Cj^kO^^a^ , in contradistinction to (the an- 
tilinear) hermitian conjugation which would have taken 
it to J, c* f.a)^d^ . It is the linearity of the transpose 
map that permits its implementation on a subsystem, the 
b-mode, of a two-mode system through 

Tr(poutdt-'d^6ti5m) ^ Tr(poutat^a'St"6')- (10) 

This partial transpose relation proves very useful. 
Beam Splitter and Conversion of Nonclassicality into En- 
tanglement: The beam splitter is a passive, energy con- 
serving system. Its unitary action U on the two-mode 
Hilbert space H"' can be viewed in two equivalent 
ways, and we will make use of both in succession. In the 
Heisenberg-type view, the effect of the beam splitter is 
to produce a 7r/4 rotation in the {a, b) plane: 

UaU-^ = ^{a + b) , UbU"^ = ^{b ~ a) , 

V 2 V 2 

U-^dU ^ ^{a-b), U-'^bU = ^{b + a) .{11) 

V 2 V 2 
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Alternatively, in the Schrodinger-type view, wherein U 
acts on the state rather than on dynamical variables, the 
fact that U commutes with the operator a^a + Wb and 
hence will not connect subspaces of the two-mode Hilbert 
space Ti*^ ® Ti.'' differing in the total number of photons 
will prove very useful |l3|. 

We will now employ the beam splitter in a simple situa- 
tion, and later extend its use to a less restricted situation. 
1. The case of antihunched input: Let us assume that 
our beam splitter is fed with a product state 



Pin = Pa® |0)fc6(0| , 



(12) 



with the b-mode in the vacuum state. After passing 
through the beam splitter we have the output state 
Pout = UpiuU^^ ■ To test the operator pout, the partial 
transpose of Pout, for positivity we use the first strategy 
(3), making the choice A = cq + Ciab . 

Arranging ci , C2 into a column C so that is the row 
vector (c*, Cj), we obtain for the left hand side of (3) 



Tr (pout^^^) = Ct 



1 



Tr(pouta^) 



C. 
(13) 

We now use the result (10) of the PT operation and eqs. 
(11), to relate these traces to corresponding expectation 
values in the original input state pa of eq. (12): 

Tr(pouta6) = Tr(pouta^^) 

= TvipinU-^b'^aU) 

= iTr(pi„(6t + at)(^_5)) 

= ^Tr(paa^a) = ^(?^a) • 

The other independent matrix element in (13) is, by (10) 
followed by (11), 

Tr (pouta^ ab^ b) = Tr (pout 0^06^6) 

= Tr{pinU-^a^b^abU) 

= iTr (pi„(at _ 6t)(at + 6t)(a - 6)(a + b)) 

= ^Tr {paO.'' a'' aa) = ^((ria^) - (?ia)) ■ 



As a result, (13) reads: 

Tr(poutAtA) = f 1/ )c. 

The determinant of this matrix is exactly the expression 
j{{Ana^) — (na))which appears in the antibunching con- 
dition (9). That is, if the input pa is antibunched, then 
Pout is necessarily NPT. We have thus proved 
Theorem T. If a product state of the form pin — Pa ® 
|0)fcb(0| , with antibunched pa, is fed into a beam splitter, 
the output is definitely NPT entangled. That is. 



((Ana^) - {ua)) < for Pa 



Pout ^ 

Pout is NPT. (14) 



2. The case of general nonclassical PND input: In the 
previous case we allowed the input a-mode state pa to 
be any antibunched state. Now we retain the product 
form (12) for the input two mode state, but choose for 
Pa a more general phase invariant state determined com- 
pletely by some given PND probabilities {p{na)}: 



Pa {{P{na)}) = V p{na)\na){na\ 



na=0 



(15) 



The probabilities p{na) are of course real nonnegative and 
normalized to unity, but are otherwise free. Given them, 
the problem of inverting eq (6) to obtain P{Ia) is one of 
the classical moment problems, in fact the Stieltjes mo- 
ment problem|lOl| as we are concerned with the variable 
la G [0,00). In this case P{Ia) is uniquely determined 
by {p{na)}] and the condition P(Ia) > translates, in 
a necessary and sufficient manner, into positivity condi- 
tions for a double hierarchy of real symmetric matrices 
of increasing dimensions, formed out of the probabili- 
ties {p{na)}. These are best described in the following 
formQ: 

Qn^ = njpina) 
Pa {{p{na)}) is a 'classical' state 

<^ the PND {p{na)} is 'classical' . 
^ P{Ia) > 

^ ^(N)^ L(^)>0 for TV = 0,1,2,... 



/ 90 qi q2 

qi <72 qs 



\ QN QN+I qN+2 

( q\ 92 q-i 
q2 qs 94 



qN \ 

qN+i 



q2N / 

9^+1 \ 

9^-1-2 

q2N+i j 



\ 9^+1 9W+2 9^+3 

Now pass the two-mode product state 

Pin — Pa ({P('^a)}) ® |0)6b(0| , 

through the beam splitter. The output state is 

Pout = UpinU^^ 

= ^(at)-|0,0)(0,0|(a)-[/-i 



(16) 



(17) 



00 



= E|£ri('^' + ^'riO,0)(0,0|(a + r 



2""nJ 

p{na)nal \r,na - r){s,na - s\ 



00 

S 2- ,:^^oV^!K-r)!.!(n,-.)! 



E 
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The general matrix element of this density matrix is 



(no ,nb \pout\na,nb) 



{ug + ni,y.p{na + Tifc) 



The matrix elements of the partial transpose pout are 
obtained simply by interchanging Ub and Ub : 



{Ua ,nb \Pout\na,nb) 



,(18) 



Now we test for possible lack of positivity of pout by 
using the second strategy (4). Consider first the principal 
submatrix of (18) obtained by setting Ub = Ua , Ub = Ua ■ 
This submatrix i?'^^ has elements 



H'^V ^ = {Ua ,na \pont\na,na) 



Clearly, this (N+1) dimensional submatrix coincides with 
modulo congruence by a diagonal matrix D^^^: 



^{N) ^ Jj{N) j^{N) j^iN) ^ 
D^-^^ ' = (2""71a!)"''J„ „ 



(19) 



Consider similarly the principal submatrix of (18) corre- 
sponding to the row (and column) choices Ub — ria +1 
and Ub = Ua + 1. This submatrix H'^^^ has elements 

H^y^^ = {na ,na + l\pout\na,na + I) 

Iria+ria' +1 



2"»+"'' +VK' + l)!rta'!(rta + l)!na! 



This (N+l) dimensional submatrix is seen to coincide 
with L^^^ modulo congruence by a diagonal matrix: 

fj(N) ^ jjiN) liN) f){N) 

= (2""+VK + l)!nj)"'<5„^,„^.. (20) 

We see that in both (19) and (20), the exhibited princi- 
pal submatrices of pout are essentially the matrices L^^^ 
and L^^) (apart from trivial congruences) appearing in 
the result (16) of the moment problem to ensure positiv- 
ity of P(/a)|3. Thus if Pai{p{na)} is nonclassical, i.e., if 
P{Ia) It Oj then one of the matrices Z^^^ will nec- 

essarily fail to be positive semidefinite, for some TV, as 
seen from (16), rendering one of the above submatrices 
of Pout nonpositive. Thus, we arrive at 
Theorem 2: If the input PND statistics {p{na)} has any 
(phase-invariant) nonclassicality, then the output two- 
mode state emerging after the beam splitter is NPT: 



Pa {{p{na)}) correspondingly 'nonclassical' 

^ corresponding submatrices of Pout t 

^ Pout is NPT entangled. (21) 
To conclude, the output nonclassicality, the transform 
by the beam splitter of the nonclassicality residing locally 
in the input a-mode photon statistics, manifests itself 
nonlocally as NPT entanglement. It is rather remark- 
able that the infinite set of matrices in (19), (20) whose 
positivity is equivalent to the positivity of pout should 
essentially coincide with the infinite set of matrices in 
(16) designed, with the help of the classical Stieltjes mo- 
ment problem 1 10l|. to capture the classicality condition 

Pila) > 0. 
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